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Abstract 

We study the finite temperature phase transition of lattice QCD with an irrelevant chiral 4- 
fermion interaction and two massless quark flavours, on 8 3 x 4 and 12 2 x 24 x 4 lattices. The strength 
of the 4-fermion interaction was reduced to half the minimum value used in previous simulations, 
to study how the nature of this phase transition depends on this additional interaction. We find 
that the transition remains first order as the 4-fermion coupling is reduced. Extending our earlier 
studies indicates that for sufficiently large 4-fermion coupling, the transition is probably second 
order. 
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I. INTRODUCTION 



In a recent series of papers [l, [2j], we have used a new staggered fermion action mod¬ 
ified to include an irrelevant 4-fermion interaction which preserves the flavour symmetries 
of the standard staggered action, to study the chiral symmetry restoring finite temperature 
phase transition of QCD. (For earlier work by others proposing similar actions see j^, [J].) 
In particular, we study the transition with two massless quark flavours. Such simulations 
are impossible with the standard action, since the Dirac operator becomes singular at zero 
quark mass. On lattices with temporal extent N t = 4, the transition appeared to be first 
order Q, while at N t = 6 it appeared to be second order as expected P|, but with the 
critical indices of a tricritical point 0, rather than the 0( 4)/ 0 (2) indices expected front 
universality arguments. We interpreted these results to indicate that at these large lattice 
spacings a, the additional interactions due to discretization, which vanish as a 2 , are large 
enough to affect the nature of the transition. For this reason, we are currently performing 
simulations at N t = 8. One such interaction is our additional 4-fermion interaction, so we 
need to vary its coefficient in an attempt to determine how important it is in determining the 
nature of the phase transition. At Nt = 6 we ran at 2 different (small) values of this coupling 
and concluded that although the j3 c = 6/g 2 of the transition changed as this coupling was 
varied, its nature did not. In this paper we present the results of studies of the dependence 
of the Ni — 4 transition on the 4-fermion coupling. 

The action we use is the staggered lattice version of the continuum Euclidean La- 
grangian density, 



( 1 ) 


For details of the staggered lattice transcription of this Lagrangian, using auxiliary fields 
to render it quadratic in the fermion fields, thus allowing simulations, we refer the reader 
to our earlier papers. The largest value of 7 = 3/A 2 used in our previous simulations at 
N t = 4, and the only value for which we determined the nature of the transition, was 7 = 10. 
We present the results of new simulations at 7 = 20 on 8 3 x 4 and 12 2 x 24 x 4 lattices, 
which indicate that the transition remains first order. We have also extended the studies of 
the 7 = 5 transition reported in our earlier work. These new studies show no sign of the 
met a-stability observed at 7 = 10, 20, suggesting that the 7 = 5 transition is second order. 
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In section 2 we present the results of these simulations. Section 3 gives our conclu¬ 


sions. 


II. SIMULATIONS AT N t = 4, N f = 2 AND 7 = 20, 5. 

We performed simulations using the hybrid molecular dynamics method with noisy 
fermions allowing us to tune to 2 flavours. We ran on 8 3 x 4 and 12 2 x 24 x 4 lattices 
to observe finite size effects. On both lattices our molecular dynamics time increment was 
chosen to be dt = 0.05, which appeared adequate. Our quark mass was set to zero and 
7 = 20. 

For the 8 3 x 4 simulations we ran for 20, 000 molecular dynamics time units for 
/ 3 = 5.27 and 7 (3 values in the range 5.285 < (3 < 5.3, and 10, 000 time units for 3 (3s 
outside this range. We observed evidence for metastability with 2-state signals both in the 
chiral condensate and in the Wilson line (Polyakov loop) for 5.285 < (3 < 5.3, with the 
clearest signals at (3 — 5.29 and (3 = 5.295. A histogram of the chiral condensate at (3 = 5.29 
is shown in figure |H showing two clearly separated peaks. From this we would conclude that 
there is a first order transition. We estimate the transition to occur at (3 — (3 C — 5.2925(50). 

For our 12 2 x 24x4 simulations, we ran for 50, 000 time units at [3 — 5.2875, (3 = 5.289 
and (3 = 5.29, close to the transition, and for shorter ‘times’ at 8 other (3 values in the range 
5.265 < (3 < 5.35. We observed clear signals for metastability with a well defined 2-state 
signal at (3 = 5.289 and (3 = 5.29, but not outside this region, clear evidence that there 
is a first order transition at (3 C = 5.289(1). The time evolution of the chiral condensate 
at (3 = 5.289 is shown in figure [21 A histogram showing the two-state signal in the chiral 
condensate at (3 — 5.289 is shown in hgureOl The separation of the peaks in this histogram is 
~ 0.56, while that on the (3 = 5.29 histogram on the 8 3 x 4 lattice is ~ 0.61. The decrease in 
going to the larger lattice is small, leading us to the conclusion that the first order transition 
is real, and not a small lattice artifact. For the Wilson line, the peak separation is ~ 0.22 
for the larger lattice and ~ 0.32 for the smaller lattice. Most of this change comes from 
a rise in the value of the Wilson line for the low temperature phase, and might well be a 
fermion screening effect on the larger lattice. 

Finally, figure EJ shows the two order parameters for the 12 2 x 24 x 4 lattice simulations 
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FIG. 1: Histogram of measurements spaced at 2 time-unit intervals on an 

8 3 x 4 lattice at /3 = 5.29. 

as functions of f3, showing how abrupt the transition is. 

Since onr earlier studies at stronger 4-fermion couplings were inadequate to determine 
the nature of the phase transition, we have extended the 7 = 5 simulations to 12 2 x 24 x 4 
lattices and added some [3 values in the crossover region to our 8 3 x 4 simulations. For 
the larger lattice we ran for 50, 000 time units for each (3 in the range 5.415 < (3 < 5.440. 
Figure 0 shows the chiral condensate and Wilson line for the 12 2 x 24 x 4 lattice. The 
transition appears smooth, and the time dependence of these order parameters shows no 
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FIG. 2: Time evolution of (V^} 2 — (tAs^sV’) 2 ) on a 12 2 x 24 x 4 lattice at /? = 5.289. 

sign of metastability for any of these (3 values. This suggests that the 7 = 5 chiral transition 
is second order and occurs at f3 = 5.425(5) 

III. CONCLUSIONS 

We have simulated lattice QCD with 2 flavours of massless staggered quarks and an 
additional chiral 4-fermion term on N t — 4 lattices. The chiral finite temperature transition 
at weaker 4-fermion coupling (7 = 20) appears to be first order, as it was at stronger 4- 
fermion coupling (7 = 10). The chiral condensate in the confined phase at the transition 
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500 


400 


m 

a 300 

CD 

> 

200 


100 


0 

0.0 0.2 0.4 0.6 0.8 

iff) 

FIG. 3: Histogram of (?/y/’) 2 — (Vh's&V’) 2 ) on a 12 2 x 24 x 4 lattice at /3 = 5.289. 



was 0.60(1) at 7 = 20, compared with 0.76(1) for 7 = 10 (in the deconhned phase this 
condensate would vanish). The change in the Wilson line across the transition is 0.23(1) for 
7 = 20, compared with 0.33(1) at 7 = 10. A priori, this softening of the transition as the 
4-fermion coupling weakens might suggest that it would become second order at vanishing 
4-fermion coupling. However, we see that at stronger coupling, 7 = 5, the transition is even 
softer and probably second order, so the situation is not so simple. 

We suggest that the order of the transition is determined primarily by the ‘standard’ 
(7 independent) part of the action, and hence by /3. On N t = 4 lattices, j3 c = 5.289(1) at 


6 




































FIG. 4: Chiral condensate and Wilson line on a 12 2 x 24 x 4 lattice as functions of (3 at 7 = 20. 

7 = 20 compared with f3 c = 5.327(2) at 7 = 10 and 5.425(5) at 7 = 5. At 7 = 10,20 the 
lower values of j3 c mean larger lattice artifacts, which we suggest make the transition first 
order. The larger /3 C at 7 = 5, has smaller lattice artifacts, which we suggest leave it second 
order. It is interesting to note that f3 c at 7 = 5 and N t = 4 is very close to /3 C at 7 = 20 for 
N t = 6 , which we have determined to be second order (tricritical), which is consistent with 
this scenario. The reason the discontinuity in the 7 = 10 chiral condensate is greater than 
that at 7 = 20 is that the stronger 4-fermion coupling enhances condensate formation. The 
larger f3 of the 7 = 10 transition increases the Wilson line just above the transition over 
the corresponding 7 = 20 value. Such a scenario suggests that at zero 
7 = 00 , the transition would remain first order. As was noted in reference 


fermion coupling 
, the presence of 
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FIG. 5: Chiral condensate and Wilson line on a 12 2 x 24 x 4 lattice as functions of (3 at 7 = 5. 

this first order transition could explain the failure of universal critical scaling at the N t — 4 
unite temperature transition for lattice QCD with 2 staggered fermion flavours as m —> 0 

MQ. 

Finally, we should comment that, because this 4-fermion term also breaks the same 
flavour symmetries as the standard staggered action, it has the potential to significantly 
increase this breaking and make improving this action more difficult. We are planning 
hadron spectroscopy calculations to determine how bad this breaking might be. However, 
this additional 0(a 2 ) symmetry breaking could be reduced to 0(a 4 ) by addition of the terms 
required to make this 4-fermion interaction that of the staggered lattice implementation of 
the SU(4) x SU( 4) or (SU( 2) x SU(2 )) x (SU( 2) x SU{2)) Gross-Neveu model 4fJ- The 
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only problem is that such an action has a complex fermion determinant, essential to the 
physics of interest. However, changing some of the couplings from imaginary to real removes 
this difficulty, and it might be possible to analytically continue results from real values to 
imaginary values of these couplings, as has been suggested for similar terms in the highly 
improved staggered actions proposed by the HPQCD collaboration [ll|J. We note that the 
4-fermion term in our action, and at least some of the additional 4-fermion terms which we 
are suggesting, are present in the HPQCD action after a Fierz transformation. 


Acknowledgments 

DKS is supported by DOE contract W-31-109-ENG-38. JBK is supported in part by 
NSF grant NSF PHY-0102409. The 7 = 20 simulations were performed on the Cray J90’s 
and Cray SVl’s at NERSC. The 7 = 5 simulations were performed on the HP Superdome 
at the University of Kentucky under an NRAC grant. 


[1] J. B. Kogut, J.-F. Lagae and D. K. Sinclair, Phys. Rev. D58, 034504 (1998). 

[2] J. B. Kogut and D. K. Sinclair, Phys. Lett. B492, 228 (2000); J. B. Kogut and D. K. Sinclair, 
Phys. Rev. D64, 034508 (2001). 

[3] K. I. Kondo, H. Mino and K. Yamawaki, Phys. Rev. D39, 2430 (1989). 

[4] R. C. Brower, Y. Shen and C.-I. Tan, Boston University preprint BUHEP-94-3 (1994); 
R. C. Brower, K. Orginos and C.-I. Tan, Nucl. Phys. B(Proc. Suppl.) 42, 42 (1995). 

[5] R. Pisarski and F. Wilczek, Phys. Rev. D29, 338 (1984). 

[6] C. Bernard et al. (MILC collaboration) Nucl. Phys. B(Proc. Suppl.) 63A-C, 400 (1998); Phys. 
Rev. D61, 054503 (2000). 

[7] G. Boyd with F. Karsch, E. Laermann and M. Oevers, Proceedings of 10th International 
Conference Problems of Quantum Field Theory , Alushta, Crimea, Ukraine (1996). 

[8] S. Aoki, et al. (JLQCD Collaboration), Phys. Rev. D57, 3910 (1998). 

[9] E. Laermann, Nucl. Phys. B(Proc. Suppl.) 63A-C, 114 (1998). 

[10] D. J. Gross and A. Neveu, Phys. Rev. D10, 3235 (1974). 

[11] Q. Mason, et al, hep-lat/0209152 (2002). 


9 






